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! Abstract 

. Rescattering effects of baryon and antibaryon in heavy quarkonium decays are investigated by 

^ ' studying their angular distributions. The rescattering amphtudes are phenomenologically evaluated 

I by modeling the intermediate range interaction as a ct or pion meson exchange between qq quarks. 

' The results show that the rescattering effects play an important role in determination of the angular 
distribution in heavy quarkonium decays. Especially, for J/ip and ip{2S) decays into AA, E"S° and 

' S^S"*" the angular distribution parameters could turn to be negative values in the limit of helicity 

\^ . conservation. These results provide us a possible explanation for understanding the negative sign 

' of the angular distribution parameter measured for J /ijj — s- YPYP , namely, it might come from the 

fH , baryonic SU(3)f symmetry breaking by incorporating rescattering effects. 

ix. 
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X ■ 1 Introduction 



Exclusive decays of heavy quarkonium into baryon antibaryon pairs (BB) are widely investigated 
as a laboratory to test perturbative QCD (pQCD) properties or to study the baryonic properties 
niElSllS]. The branching fractions of J/^p and ip' decays into BB have been systematically studied 
in the framework of pQCD The angular distributions of e+e~ J/tp, tp' B%]3% {B%: octet 

baryon) , were also investigated theoretically and experimentally by many groups |S1 IHl El E| ■ 
For J /ip produced from e^e~ annihilation, they must have spin ±1 along the beam direction. For 
exclusive decays of J /ip into octet baryonic pair, the outgoing BB can be taken either S-wave or D- 
wave state in laboratory system, thus the experimental observation of the outgoing baryon direction 
takes the form dN/d cos 9 cx (1 + acos^ 9), where 9 is the angle between outgoing baryon direction 
and e+e~ beam direction. The angular distribution parameters of a have been measured by many 
groups |1UI nil [T^ . Recently, the BES collaboration have released a more accurate measurement 
on the parameter a for J/i/i —> pp, AA and E'^E" decays with 58 million J/tp events accumulated 
at BEPC as shown in Table It is interesting to note that the parameter a for J/ip ^ E^E" is 
negative within la accuracy, and the sign of the central value of a for J/ip — > measured at 

Mark II is also negative though its uncertainty is still very large. 

Theoretically, the study on the angular distribution parameter provides us a copious information 
on polarization of decaying particles , decay mechanisms, and strong interaction and so on. Brodsky 
and Lepage once demonstrated that the exclusive process at large momentum transfer can be 
used to test the gluon spin and other basic elements of perturbative quantum chromodynamics 
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Tab. 1: Angular distribution parameter a for J/t/; BB decays. They 
are assumed to be the form of dN/d cos 9 cc 1 + a cos^ 6. 



Calculated value of a 



Decay mode 


Measured value of a 


Ref. p 


Ref. m 


J/ip pp 


0.68 ± 0.06 j9j 


1 


0.69 


J/ijj aa 


0.65±0.1l[in] 


1 


0.51 




-0.24 ± 0.20[in] 


1 


0.43 




-0.13±0.59[Il] 


1 


0.27 


ip' — > pp 


0.67 ± 0.16(12] 


1 


0.80 



(QCD) They verified that the vector-gluon coupling conserves quark helicity when quark and 
gluon masses are negligible at the large scale . Since the hadronic helicity is the sum of the 
helicities of its valence quarks, thus leads to the conclusion that the total hadronic helicity is 
conserved up to corrections of order viq/Q or higher. For e+e~ — > J/ij) — s- BB decays, angular 
distributions take the form of oc 1 + cos^ 9 by the selection of the helicity conservation rule, 

which gives the asymptotic value of a — 1. They conclude that verifying the 1 + cos^ 9 angular 
distribution would provide one with a clear proof of the spin one assumption for the gluon, since 
with a scalar or a tensor gluon one would be led to a sin^ 6 distribution. 

In the past, this asymptotic value has been corrected by many groups. In charmonium decays, 
we are far from an asymptotic regime due to the factor that tob/to^ ^ 1/3 is far from negligible. 
The baryonic mass effects have already been investigated by Claudson, Glashow and Wise jj], and 
they found a value of a = 0.46 for J/ip ^ pp, which confirms the idea of large mass corrections 
to the asymptotic value. The quark mass effect has been investigated by Carimalo It is 
assumed that the quark created with effective mass to inside a baryon carries one-third of the total 
momentum of the out-going baryon, and the three symmetric gluons carry the same momentum 
(M^/3) with no other extra interaction. In this simple static-quark model, the angular distribution 
can be explicitly expressed by 



a — - -T^ ; with u= — f-. (1) 

(1 -t- uY + u(l + 6-it) to| 

This formula gives the parameter a — 0.69 for J/^ —f pp, very close to the experimental value. 
However, the calculated branching fraction of J —> BB decays are inconsistent with measured 
values with Icr in this model. Some authors argued that the non-perturbative effects of the baryon 
should be parameterized in the calculation. Further investigations including quark mass effects, 
electromagnetic effects and higher-twist corrections to a have been made with a general conclusion 
of < a < 1 |S1|H|. 

Phenomenologically, some authors argued that, in charmonium decays, rescattering effects be- 
tween the final particle states sometimes are essential to explain experimental observations. For 
example, a pp rescaterring near their threshold could make some contributions to the observation 
on the near threshold-enhancement of J/ip —> jpp ^ . In meson sector, one example is that the 
rescattering of a2p and aip into pir can change the pn production rate substantially, which could 
be one of explanation for the pir puzzle [Jj. Whether the BB rescaterring could change their dis- 
tribution significantly? Stimulated by recent BES observations, we make a further investigation on 
the angular distribution for heavy quarkonium decays into BB by taking rescattering effects into 
account. 



2 Model and Formulation 

The BB produced from heavy quarkonium decays moves oppositely with higher momentum (see 
Figure n (a)), they might experience an intermediate-range interaction by exchange of mesons 
(Figure n(b)). We assume that the outgoing baryonic pair survives from the interaction without 
annihilation. We only consider one-pion and tow-pion exchange processes. In principle, other 



mesons with higher masses might be exchanged, their contributions are supposed to be less than 
pion-exchange processes. As traditional treatments to the NN interaction, the process of two- 
pion exchange is described by the picture of a meson exchange |15|. For consideration of isospin 
conservation, one-pion exchange could be allowed in pp and S^S"*" decays suppressed by a isospin 
factor of 1/3, and not allowed in AA and S^E^ decays. As in chiral quark model we model the 
intermediate range attraction in the BB interaction as a rescattering of qq quarks by exchanging 
a (T or pion meson (Figure ^(c)). 




Fig. 1: The schematic illustration for heavy quarkonium decays into BB (a) and 
via a rescattering process of BB (b), which is equivalent to a exchange of a meson 
between qq quarks (c). 



The helicity amplitude corresponding to Figure ^ (a) and (b) can be expressed as: 

Ml = ^lUi-Ai {d2,h)Fx'„y, , for FigureH (a) (2) 
M2= E (2J+l)^lU-A. (^i>i)^A,,A.r/.,-, /?3::-A..A;-Ai (^2, <A2), for Figure □ (b) (3) 

Ai,A2,J 

where M,Xi{X[) and A2(A2) denote the helicity of heavy quarkonium, baryon and antibaryon, 
respectively. For heavy quarkonium production from positron annihilation, M = ±1. F and T 
are helicity amplitudes corresponding to heavy quarkonium decays into BB and BB rescattering 
processes, respectively. For considering parity invariance and time reversal invariance, one has 

F — = = F+_,T_+ — = -T + = -r+_++ = T+++_ T+ -T — ++ = 

-T_+++ = T++_+ = ri,T++__ = = r2,r+_+_ = = t_++_ = t+__+ = and 
T++++ = T___ = n. 

Note that interference terms vanishes between Mi and M2 due to the orthogonal D— function, 
then the angular distribution of the decayed baryon is expressed by: 

dV 



dcos02 



J2 \Mi\^d(j)2+ \M2\^deid(j)id(j)2 0il + a'cos^e2, (4) 



with 



a' = {4a - 97r2[(a - 1)(4T2T4 + 2T|) + 6Tf (a + 1) + a{2Tl - 3T|) - 2T| - 3r|]} 

/[^ + 27^\2Tl + Ti){l + a)], (5) 

where a = ^^^U , „ L'^^U is the net angular distribution parameter for BB without rescattering 



effects. If the rescattering effects vanish, i.e. Ti —i- 0, then a' = a. It is worthwhile to note that 
the sign of a' could be positive or negative determined by BB scattering amplitudes. 

Unfortunately, evaluation of rescattering amplitudes in the framework of pQCD cannot be made 
without models, this is mainly due to the fact that the non-perturbative information on bound states 
can not be fully determined from a reliable method in theory. For low lying states of octet baryons, 
we will employ the so-called constituent quark model to evaluate rescattering amplitudes. For a 
exchange, we define an operator 

xS'^{pi — p'i)S'^{p2 — P2)^'^iP3 ^ P3 ^ 2;) + permutation terms]. 



where is the effective couphng strength of a a meson to qq quarks, which will be fixed by using 
the angular distribution parameter of J/ip ^ pp as input value. and mo- are the cr's decay width 
and mass, respectively, x denotes the transferred momentum. u(j)i,Si) and v{qi,Si) are Dirac 
spinors for quarks and antiquarks with helicity Si and Si, respectively. Their normalization is taken 
as uu = —vv = m/E. For one-pion exchange, the should be replaced by 75 (/tt and propagator 
by l/(x^ — + ie). Then the rescattering helicity amplitude is evaluated by: 

1=1,3 

+^AiA2,a;a^ + -4aiA2,a;,a; 1-05(^1, Si)) 

xS^P -Pi-P2~ P3)6^P' -1^,-1/2- I^3)^'{Q - qi - 92 - q3)S'{Q' - (f, - - q^i). (7) 

where ipsiPi, Si)ip§{qi,Si) is the product of B and B wave function in momentum space, which 
includes the spin, flavor and spatial wave function. They are constructed in the naive quark model. 



3 Numerical Results 

3.1 SU(6) basis 

For illuminating the role of the effect of SU{3)f symmetry breaking in rescattering processes, we 
firstly evaluate rescattering amplitudes by using SU(6) basis and then compare it with the so-called 
uds basis. We naively assume that the masses of quarks satisfy m„ = rrid — rng = ms/S, where 
niB is the baryonic mass. This means that we neglect binding energies in the baryonic bound state 
and use the non-relativistic static quark model. The spin-flavor wave functions of octet baryon are 
described by: 

0|^ = 0|j. = -i=(x''0'' + xV), (8) 
where x'' and are the mixed-symmetry pair spin-i wave function. For example, we have 

x?/24/2 = -^{inT> + im>-2im>}, 

Xi/2,1/2 = ^{inT>-UTT>}, (9) 

for the case of the total spin i and its projection i. The flavor wave function and 0'^ are exactly 
analogous to the spin wave function but in the flavor space. As in the most naive quark model, 
we assume that the bound state wave function for constituent quarks can be described by a simple 
harmonic-oscillator eigenfunction in their center-of-mass (cm.) system, i.e. 

0(fcp,fcA)=(^e-(^^^^l)/^^ (10) 

where (3 is the harmonic-oscillator parameter and kp, k\ are defined as kp = {ki + k2 — 2k3)/^/6, 
and kx = (fci — k'i)l\f2 with fci, ^2 and fca the momenta for the three quarks in cm. system of the 
corresponding baryon. 

3.2 uds basis 

To count for the SU(3)f symmetry breaking in strange baryons, we use the so-called uds basis, 
which makes explicit SU(3)i? symmetry breaking under exchange of unequal mass quarks [l7| . 
Flavor functions for strange baryons are taken as 

0A = —j=(ud~du)s, (11) 
v2 

0s ^ -^{ud + du)s, (12) 

02 — ssd, (13) 



and the construction of the spin wave function x proceeds exactly analogously to that of the flavor 
wave functions. For instance, the spin wave functions for the total angular momentum quantum 
1/2 with their z projection 1/2 are: 



XA = ^(inT>-im», 

xs = ^(1 TiT> + UTT>-2| TTi», 
XH = ^(ITiT> + UTT>-2|TTi». 



(14) 
(15) 
(16) 



The spatial wave function for the ground state of the strange baryons are chosen similarly to those 
of protons, but have an asymmetry between the p and A oscillators. In the baryonic cm. system, 
it reads: 

1 



(^2/3p/3^)3/4 



,-(fe^/2/3p+fel/2/3;,) 



(17) 



where Pp = (Sfcm)^/^ and m is the quark mass. For example, we take m = m„ for A, S, and 
take m = rris for S. (3x = (Sfcm)^/^ with m\ = 3mTO3/(2m + 1713) > m {m^ = nis for A, S, and 
7713 = m„ for S). Using these equations, we relate /3p(A, E,S) and /3a(A, E,S) to the harmonic- 
oscillator parameter (3 of nucleons, i.e., /3p(A, S) = (3, (3p{^) ~ y/ms/md(3, and /3a(A, = 
y/mx/mPp{A,T,,E). 

Tab. 2: The angular distribution parameters for J/ip, ip' and T(IS') decays into pp, AA, 
S^S" and incorporated BB rescattering effects by exchange of a o" meson. The 

parameters are set as m„ = rrid = = ms/'i and (3 = O.lGGeV^ for SU(6) basis, while 
mu = md = 310MeV,ms = 490MeV and /3 = O.lGGeV^ for uds basis. 



Channel 



J/^p pp 
AA 



SU(6) basis 

y _ 4a+(14.8-8.5a)g^ 

4+4.6{l+a)g4 
, ^ 4a+(14.8-8.3a)g4 
4+4.9{l+a)g4 
4a+(13.5-7.6a)3^ 

4+4.6{l+a)g4 
4a+(8.7-4.9a)g4 



uds basis 



4a+(14.8-8.5a)g^ 



a 



a 



ip' pp 
AA 



a 



a 



a 



T(15) pp 
AA 



a 



a 



a 



a 



4a+(14.8-8.9a)g^ 
■ 4+4.3{l+a)g4 
4a+(18.0-10.5Q!)g^ 

4+5.5(l+Q)g4 
4a+(19.4-ll.la)g^ 

4+6.0{l+a)g4 
4a+(20.4-11.4a)g^ 
4+6.6(l+a)g;^ 
_ 4a+(1.4-1.0a)g| 
" 4+0.3{l+Q)g^ 
_ 4a+(2.1-1.5a)g^ 
" 4+0.39(l+a)g4 
_ 4a+{2A-1.7a)gi 
- 4+0.5{l+Q)g4 
_ 4a+(3.1-2.1a)g4 
" 4+0.6(l+a)g4 



4+4.6(l+a)g4 
/ _ 4a+(13.1-164.9a)g^ 
" ~ 4+118.1(l+a)g4 

, _ 4a+(11.6-25.7a)g^ 
" 4+20.8(l+a)g4 

/ _ 4a+(7.0-15.4a)g4 
" - 4+12.5(l+a)^ 



a' 



4a+(14.8-8.9a)g;^ 



a 



4+4.3{l+«)gJ 
, _ 4a+(176.6-207.9a)g^ 



4+139. 
, ^ 4a+(18.6-36.6a)g4 

4+30.2(l+a)g4 
, ^ 4a+(15.7-32.3a)g4 
" 4+26.5(l+a)g:^ 



y _ 4a+(1.4-1.0a)g^ 

4+0.3(l+a)g4 
, ^ 4a+(25.6-26.8a)g4 
4+9.9{l+a)g4 
4a+(2.7-4.5a)g^ 
4+3.5(l+a)g4 
/ _ 4a+(1.7-3.4a)g^ 
4+2.6(l+a)g;| 



a 



a 



For numerical evaluation of rescattering amplitude of BB^ there are three parameters to be 
determined, i.e. the quark mass rriu, and the harmonic oscillator parameter f3 in nucleon 
spatial wave functions. As usually used in naive quark model, we set m„ = uis = 310 MeV, 
rUs = 490 MeV and /3 = 0.16 GeV^. The rescattering amplitudes Ti{i = 1,2,3,4) are obtained by 
evaluating the integration of Equation {Tj) directly. The numerical program can be obtained from 
authors. Then we list the angular distribution by incorporation of BB rescattering effects in Table 

m 



The observed angular distribution of BB depends on the effective coupling strength g^, and the 
net parameter a. In the chiral quark model, the coupling strength g is determined from gT^NN 
which is g — 2.62. Here we fix this parameter by using the experimentally angular distribution 
of J jt^} —f pp as an input value. For consistency the net parameter a should be calculated from 
the same quark model. Carimalo formula is a good approximation in the description of the net 
parameter a by incorporating baryonic structure information. If the quark mass is negligible or 
4to^ ^ (M : heavy quarkonium mass), the net parameter tends to the asymptotic value 
a = 1. We will see how the observed angular distribution parameter changes when the quark mass 
vanishes, so as to check its asymptotic behavior under the condition of the helicity conservation. 

As shown in Tabled if the difference between the light quark masses and strange quark masses is 
ignored and octet baryons are described by SU(6) wave function, the observed angular distribution 
parameters remain positive when the net parameter a tends to the asymptotic value. Though 
the experimental value for J/ip AA seems to be consistent with results from SU(6) basis, it is 
not enough to distinguish the SU(6) basis from the uds basis due to the uncertainty of the net 
parameter a. It is worthwhile to note that calculated values using the uds basis for J/ip decays 
into strange baryons change their signs when the a tends to the asymptotic value. Though the 
experimental error for J/ip — > S''li° is still larger, the measured central value for this decay is 
negative within la accuracy. So this result provides us a possible explanation for understanding 
the negative sign of angular distribution parameter for J/ip TPT,^ decay, namely, it might come 
from the baryonic SU(3)f symmetry breaking by incorporating rescattering effects. For further 
test of rescattering effects in charmonium decays, higher accurate measurements of the angular 
distribution for J/tp,ip' S^E, are desirable. 



Tab. 3: Comparison of the angular distribution parameter between experimental and theoretical 
values by SU(6) basis and uds basis. Where the values of a'{a = 1) and a' (Carimalo) indicate 
that the net parameter a are set to one and evaluated by Carimalo formula, respectively. 





SU(6) basis 


uds basis 


Experimental value 


Channel 


a'{a = 1) 


a'(Carimalo) 


a'{a = 1) 


a' (Carimalo) 


a' 


J/^jJ pp 


0.68(input) 


0.68(input) 


0.68(input) 


0.68(input) 


0.68 ± 0.06 ^ 

0.65 ±0.11 [ini 


AA 


0.66 


0.69 


-0.13 


0.37 




0.65 


0.63 


-0.34 


0.23 


-0.24 ±0.20 Uni 




0.64 


0.45 


-0.33 


0.18 


-0.13 ±0.59 IT] 


ip' — > pp 


0.70 


0.84 


0.70 


0.84 


0.67 ±0.16 Jij 


AA 


0.69 


0.80 


-0.11 


0.20 






0.68 


0.78 


-0.30 


0.19 






0.68 


0.75 


-0.31 


-0.03 




T{lS)^pp 


0.84 


0.98 


0.84 


0.98 




AA 


0.80 


0.97 


-0.06 


0.59 






0.80 


0.96 


-0.24 


0.74 






0.78 


0.95 


-0.31 


-0.28 





The case of one-pion exchange is also investigated in pp and S~S+ decays under the SU(6) basis 
and uds basis. Parameters are chosen as same as the cr-exchange case. Numerical results show that 
the contribution of one-pion exchange from the two mixed-symmetry spin wave functions plays a 
different role in the evaluation of rescattering amplitudes. Since the coupling of qqir is selected 
8^ ffTrTS) the contribution from the A— type spin wave function is dominant over type functions. 
Thus the values of rescattering amplitudes for S~S+ decays under the uds basis are quite smaller 
than those under the SU(6) basis. With the same treatment as to the a exchange case, the effective 
coupling strength of g^^ is determined by using the experimental value of J/^ — > pp decay as an 
input value, we find out an approximate relation a' ~ a for J/V', V'(2»S') and T(IS') decays into pp 
and channels under the two bases used. 



4 Discussion and Summary 



Based on our evaluation, we suggest that charmonium decays into I]°fi° and S^S+ could be used 
to look for BB rescattcring effects. For T(IS') — > BB decays, the branching fraction is relatively 
smaller than that for charmonium decays, furthermore, its BB rescattering amplitudes are largely 
suppressed by a factor of l/Af^. It should be pointed out that the BB angular distribution dose 
not become isotropic at the limit of the baryonic zero- velocity, instead it is dominantly determined 
by rescattering processes. From the point of this view, the rescattering effects should be looked for 
in the decay with a slowly moving BB pair. 

It is instructive to compare our results with ones from the quark-scalar-diquark model in |18| . In 
the asymptotic limit, i.e., the helicity conservation, the sign of parameter a' from uds basis turns to 
be negative. In |18| the negative sign of this parameter appears when the masses of quarks, gluons 
and baryons tend to zero for point-like baryons. This implies that if some interactions between 
out-going quarks are considered, the sign of the angular distribution parameter could flip in the 
asymptotic limit of the helicity conservation. 

To summarize; the rescattering effects of BB in heavy quarkonium decays are investigated by 
studying the angular distributions. The rescattering amplitudes are phenomenologically evaluated 
by modeling the intermediate range interaction as a cr or pion meson exchange between qq quarks. 
The results show that the rescattering effects play an important role in determination of the angular 
distribution in heavy quarkonium decays . To compare the calculated values by using the SU(6) 
basis and the uds basis, one finds that the parameter a' turns to be negative values by suing 
the uds basis in the limit of asymptotic value a = 1. Even though the accurate values of a' are 
not fully determined in our calculation due to uncertainties of the net parameter a, we conclude 
that these results still provide us a possible explanation for understanding the negative sign of 
angular distribution parameter measured in the decay J/?/; — > E^S", namely, it might come from 
the baryonic SU(3)p symmetry breaking by incorporating rescattering effects. For further test of 
rescattering effects in charmonium decays, higher accurate measurements of the angular distribution 
for J/'0, 4'' ~^ S^S, 5~S+ are desirable. 
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